Microscopic theory of the Andreev gap 
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We present a microscopic theory of the Andreev gap, i.e. the phenomenon that the density of 
states (DoS) of normal chaotic cavities attached to superconductors displays a hard gap centered 
around the Fermi energy. Our approach is based on a solution of the quantum Eilenberger equation 
in the regime ir> <C £e, where to and tv are the classical dwell time and Ehrenfest-time, respectively. 
We show how quantum fluctuations eradicate the DoS at low energies and compute the profile of 
the gap to leading order in the parameter to/tE- 

PACS numbers: 03.65.Sq, 03.65.Yz, 05.45.Mt 



The attachment of a superconductor to a conducting 
cavity leads to a suppression of the normal density of 
states - the proximity effect. For cavities with classically 
chaotic dynamics, a discrepancy is found between semi- 
classical calculations [T] and such based on random ma- 
trix theory (RMT) [2j: Semiclassics obtains a small yet 
finite DoS for all excitation energies e above the Fermi 
level cf, while RMT predicts the formation of a hard gap 
below some energy e* . The origin of this so-called 'gap 
problem' in Andreev billiards was pointed out by Lod- 
der and Nazarov some time ago pQ : quantum corrections 
not captured in the principal semiclassical approxima- 
tion are expected to generate a hard spectral gap for 
trajectories longer than the Ehrenfest time. Although 
various semi-phcnomcnological realizations of this mech- 
anism have been formulated, a fully microscopic theory of 
gap formation is outstanding. The construction of such 
a theory is the goal of the present paper. 

Quasiclassical Eilenberger equation — Consider a two- 
dimensional Andreev billiard, i.e. a chaotic normal- 
conducting cavity attached to a bulk superconductor. 
We wish to compute the cavity DoS in a 'semiclassical' 
regime where the quantum time scales of the problem ex- 
ceed all classical scales. Under these circumstances one 
expects [1] the gap, e* to be set by the inverse of the 
Ehrenfest time, e* = nh/2t E , where f E = A -1 ln(c 2 /h), A 
is the dominant Lyapunov exponent of the system, and 
c 2 a classical action scale whose detailed value is of little 
relevance. Heuristically, tg is the time a minimal wave 
package needs to spread over classical portions of phase 
space; the dynamics at time scales beyond is no longer 
classical. 

To compute the DoS, we start out from the quantum 
Eilenberger equation (for notational convenience we sup- 
press the infinitesimal imaginary increment in e + iO) 



[ea 3 - iA<7 2 + H1*G]=0 



(1) 



for the quasiclassical retarded matrix Green function, 
G(x), i.e. the Wigner transform of the Gorkov super- 
conductor Green function. In (JT|), o~i are Pauli matri- 
ces acting in particle-hole space, x = (q, p) T is a phase 



space point in the shell of constant energy, H(x.) = ep, 
-ff(x) is the Hamilton function, and the order param- 
eter amplitude A = A(q) is non-vanishing only at 
the cavity-superconductor interface. The Green func- 
tion is subject to the nonlinear constraint G * G = 1, 
and yields the DoS as v(e) — ^ Re J d 2 xtr [G(x)cr 3 ], 
where il = j H ^- ep d 2 x 1 is the volume of the energy 
shell and vq the normal metallic DoS. Finally, the sym- 
bol '*' indicates that all products between phase space 
functions in Eq. ([!]) are Moyal products (A * B)(x) = 
exp(f9£/0 x )| x=x/ A(x')B(x). 




FIG. 1: Inset: classical trajectory connecting the supercon- 
ductor/normal conductor interface of a chaotic Andreev bil- 
liard with itself. Main part: abstract phase space representa- 
tion of that trajectory and its vicinity in a system of locally 
stable (s) and unstable (u) coordinates. The meaning of the 
shaded areas is explained in the main text. 

Classical evolution and its inconsistency — Upon Tay- 
lor expansion to lowest orders (A * B)(x) = A(x)B(x) + 
f{A,B}{x) + 0{h 2 d 2 x ) ({,} is the Poisson bracket) 
Eq. Q assumes the standard form of the classical Eilen- 
berger equation [3] 



[iea 3 + Aa 2 ,G}-hCG = 0, 



(2) 



where C = {H, .} generates the classical Liouville flow. 
However, (finite order) Taylor expansions of the Moyal 
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product become problematic in cases where the function 
G displays structure on linear scales < 0(h) and higher 
order derivatives 0(h 2 d 2 ) become of the same order as 
HC; as we shall see, this is precisely what happens on 
the solutions G supporting the DoS in the region of the 
spectral gap. 

The classical Eilenberger equation ([2| describes the 
evolution of G along individual classical trajectories 7 be- 
ginning and ending at the superconductor interface (cf. 
inset of Fig. [I]) Parameterizing a trajectory 7 of length 
T in terms of a coordinate t S [— T/2, T/2], the Liouville 
operator on 7 assumes the form C = d t and the solution 
in the asymptotic limit e/A — > is [1J 
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(3) 



Denoting the <x;-components of G by Gi 
obeys the boundary conditions [T] 



the solution 



Gr,i(±T/2) = ±»G Tl 3(±T/2), G T , 2 (±T/2) = 1. (4) 

The component Gt,3 = const, generates (via the identity 
Im (tan (a; + iO)) = ^J2m S(x — (m + 1/2) it)) a quanti- 
zation condition, eT = (to + ^)ttH, m = 0,1,2,..., for 
the flight times of trajectories contributing to the DoS at 
energy e. The exponential sparsity of trajectories with 
T » <d much larger than the average dwell time H] 
then leads to an exponential suppression of the DoS for 
e < Ht^ 1 , but not to a gap. 

In view of the continuity conditions underlying the ap- 
proximation ([2]), it is mandatory to explore what happens 
as we transversally depart from an isolated trajectory 
into surrounding phase space. To this end, it is useful 
to interpret each trajectory as element of a corridor or 
band [6j [Tj [S] which is formed by all trajectories that 
run through the same sequence of scattering events. A 
schematic of a band is shown in the bottom part of Fig. [T] 
where the straight line represents a trajectory beginning 
and ending at points x and x' in the SN interface. We 
introduce Poincare sections through the trajectory, and 
span them by the locally stable and unstable coordinates, 
s and u, respectively. The shaded areas then represent 
the SN interface (Si), the image of that area under the 
Hamiltonian flow after time t (S2), the intersection of the 
image with the interface (S3), and the pre-image of the 
intersection (S), respectively. Points in S remain com- 
pactly confined and exit at the same instance T. The 
image of S under evolution defines a 'corridor' of sections 
across which the quasiclassical solutions Gt is nearly con- 
stant. While the transverse area, us, of the corridor is 
a conserved quantity, its shape is not. At a given in- 
stance of time, t, its smallest linear extensions is given 
by (cf. Fig. |l( ~ const, x min( e - A ( T / 2 +*\ e -^ T ^), 
with a classical proportionality constant. For trajectory 
times T > t^, that scale may shrink below 0(h), and 



this is when Eq. ([2| becomes problematic: at low ener- 
gies, e ~ feg 1 , the narrow corridors of long trajectories 
T > ts meander through the bulk of phase space, in 
which trajectories are of average length ~ trj <C £e an d 
Green functions are 'locked' to the superconductor or- 
der parameter, G(e) ~ a%. (Here and throughout, we 
use the notation ~ to indicate equality up to inconse- 
quential corrections scaling with some positive power of 
h.) The ensuing sharp variation of the solution Gt over 
trans-corridor sections of quantum extension < 0(h) con- 
flicts with quasiclassical smoothness conditions required 
for Eq. (|2j. 

Our solution to the problem proceeds in two steps: we 
first transversally extend ^ to a solution of ([2| in a 
'Planck tube' [5] 

Z=[jZ t , Z t = {(u,s,t)\ \us\<h, |u|,|a| <c}, (5) 
-¥<*<¥ 

centered around 7. This - singular - configuration will 
then be the basis for the construction of a smooth con- 
figuration G that solves the quantum equation up to 
corrections ~ tp /tg . The quantum G displays a hard 
spectral gap. 

Consider, then, the corridor carried by a trajectory 7 
of length T > tg. (In the wide corridors of shorter trajec- 
tories the Green function does not depend noticeably on 
transverse coordinates (u,s,t) £ Z t and the solution ^ 
can be taken face value.) We assume the corridor sections 
Zt to be small enough to afford a linearization [10 



C = d t + Xud u — Xsd s 



(6) 



where the terms ud u and sd s describe the divergence and 
contraction of phase flow around 7, respectively. 

Going forward (backward) in time, the trajectory 
through a point (u, s, t) <E Z t will stay in the vicinity 




FIG. 2: Vicinity of a long trajectory in a system of locally sta- 
ble (s) and unstable (u) coordinates. Points at the boundary 
I its I ~ h belong to trajectories 7 generically of length T ~ £e 
(here illustrated by the curved line.) The cloudy region at the 
ends of 7 represent generic phase-space points. 
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FIG. 3: On the length of trajectories piercing the boundary of 
the Planck cell around a long reference trajectory, a) reference 
times corresponding to a bulk point in the N-region, b) point 
close to the exit into the superconductor, c) point close to the 
entrance into the N-region. Discussion, see text. 

of 7 for a time t(u) (t(s)) where t(x) = A -1 ln(c/|cc|) . (cf. 
Fig. [2]) Thereafter a classically short time, typically of 
0(tr>), passes before the departing trajectory exits; up to 
classical corrections, the time of flight of the trajectory 
through (u,s,t) thus reads t(u) +t(s) = A -1 ln(c 2 /|us|). 
Specifically, for phase-space points on the boundary of 
the Planck cell \us\ ~ H and therefore t(u) + t(s) ~ ig. 
The above consideration applies to phase space points 
far away from the SN interface (cf. Fig. [3] a)). For 
points close to the interface, it may happen that the 
trajectory through (u,s,t) hits the interface before it 
has diverged up to c, in which case the exit time is 
shorter than t(u) (Fig. [3] b)). Or, it has been in the 
system for a time shorter than t(s) before the reference 
point is reached (Fig. [3] c)). We subsume these dif- 
ferent cases, by introducing effective in- and out-times 
t a (u) =min(t(u),T/2-t) and U(s) = min(*(s),T/2+t), 
where the function min(i, t') = — ^ln(e _At + e~ A * j 
smoothly interpolates between t and t' over a 'micro- 
scopic' switching interval ~ X . These functions evolve 
uniformly, in the sense £(ij/ (s/tt)) = (+/— )1. This 
means that the effective (up to corrections of 0{tr>/t-Eji) 
duration of the trajectory through (u, s, t) G Zt, is given 
by Tt u _ s t > = t (u) + ti(s) and the trajectory parameter 
by T (u,s,t) = - t (u)). Substitution, T -> T(„ jS , t ) 

and t — > T( n _ s t ) in (|3j) then obtains a transverse extension 

G c ( U , S ,i) = G T( „, s , t) (r KM) ) (7) 

of ([2]) [12 . G c solves the Eilenberger equation in direct 
consequence of the flow-uniformity of £j(s) and t (u). By 
the same token, however, the solution becomes singular 
at times T > t E when t j>0 begin to display structure on 
scales < h. Next, we show that this is not what happens 
in the full quantum dynamics. 



Quantum evolution and spectral gap — Let us define a 
generalization, tj(u,s,t), of U(s,t) by requiring unifor- 
mity under the full dynamics, — ihr 1 [H* tj] = 1, or 

£i? + [V?*?]=l, (8) 

where [V* ] = — ihr 1 ^* ] — C accounts for quantum cor- 
rections to the linearized classical dynamics. The above 
equation may be solved by introducing 'action-angle co- 
ordinates' I = us, 4> = |ln(u/s) in terms of which 
£ = dt + \d,j,. The V-term may now be formally removed 
by 'gauging' Eq. (JsJ with 

U{I^t)=Vel^ d ^\ (9) 

where el" is defined by a Moyal series expansion in the 
exponent and V is a (^-ordering prescription (see Ref. 
for details) accounting for the non-commutativity of 
V(I,<p) at different values of (j>. By construction [13], 
U obeys CU = V*U, which means that Eq. ^ is solved 

by 

tf(u, s, t) = (U*U* U~ l ){u, s, t). (10) 

In practice, both the detailed form of V and of U will 
not be known. This lack of knowledge, however, is not 
of essential concern to us; to the logarithmic accuracy 
required by the present analysis, basic scaling arguments 
suffice to determine the action of U on f j : describing non- 
linear corrections to the flow, the expansion of V for small 
u, s starts as HV = us x 0(u n s m ),n + m > 0. Accord- 
ingly, U = l+h~ 1 usxO(u n s m ). This entails that for any 
function / that is smooth (analytic) around u = s = 0, 
(Cr*/*^- 1 )^,*) = f(u,s)+0(d u fu n+1 s m ,d s fu n s m+1 ). 
At the small values of coordinates we are interested in, 
\us\ ~ h, the 0{. . . )-terms become irrelevant, which re- 
flects the irrelevancy of dynamical corrections to the lin- 
earized flow close to the trajectory center. To explore the 
effect of U on singular functions (such as U(s) ~ ln(|s|)), 
we notice that for arbitrary /(s), 

e- iku *f(s) = f(s + hk) *e~ lku . (11) 

This identity suggests to introduce a Fourier mode de- 
composition U(u, s, t) — J dkU( St t){k)e~ lku . Specifically, 
let us consider values |s| ~ h, where singularities be- 
gin to put the semiclassical theory at risk. For these 
values, the support of the mode coefficients U( s ^(k) ex- 
tends up to 'classical' values k ~ h° [TJ. We thus ob- 
tain tf(u,s,t) — J dkti(s + hk) * F( a ,t){k)i where the 
positive indefinite but normalized (J dkF( s t )(k) = 1) 
'weight' function F s (k) = {e~ lku U { ^ t) {k)) * IJ-^s,^. 
A straightforward estimate now shows that for asymp- 
totically small h the integral evaluates to tj(u,s,t) — 
ti(\s\+H/u ) — min(ti(s), t^), where u is anon-universal 
constant. Similarly, t^(u, s, t) ~ min(i (u), <e)- Summa- 
rizing, we have found that the operators U act to truncate 



singularities in trajectory times in a manner independent 
of the detailed form of the potential V. 

Building on these results it is now straightforward 
to construct a smooth solution of the quantum equa- 
tion Q: its general solution is given by G03 + (1 — 
G 2 ) (cos a 2 +sin where r q = (if - tg)/2 

and we used that smooth functions /(r 9 ) ('sin', 'cos', 
etc.) evolve linearly, [F* /(r«)] ~ f(T q )[H* T q ] = 
ihf'{T q ), up to corrections of 0(l/t E A) 15]. The nor- 
malization function C — C(u,s,t) is determined by re- 
quiring stationarity [H* C] = 0, and compatibility with 
the boundary conditions Q. These two conditions lead 
to the identification C = it&n(eT q (u,s,t)/h), where 
T q (u,s,t) = min(T(u, s, t), t&) is the effective trajectory 
time. In conclusion, we have found that the quantum 
equation (fll is solved by GTi( u .s,t){T q { u , s ,t)), which dif- 
fers from the solution of the classical equation Q by an 
upper cutoff t E limiting both the trajectory time T q and 
the trajectory parameter r q . Technically, this is the main 
result of the present letter. 

The above solution signals that quantum fluctuations 
couple narrow bands of transverse extension < h to 
neighboring phase space. This coupling is strongest 
in the terminal regions of long trajectories |r| > t E 
where bands flatten in one direction. Inspection of 
Fig. [2] shows that the ^-neighborhood of these segments 
is pierced by trajectories whose length and parameter 
are uniformly given by T ~ t E and |r| ~ f E /2, respec- 
tively. At these values the solutions G are nearly sta- 
tionary (up to corrections 0(£rj/i E ), and this reflects 
in the asymptotic constancy of the regularized solution 
G T q> tE (T q > ±i E ) — G tE (±i E ) at large parameter val- 
ues. The capping of trajectory times T q < i E in turn 
implies a vanishing of the DoS for energies e < e*. (Tech- 
nically, Re(G3(x)) = — Imtan(eT l? (u, s,t)/H) is vanishing 
for these energies.) The fact that all trajectories of nom- 
inal length T > tft get reduced to the uniform effec- 
tive length T q ~ i E implies an accumulation of spectral 
weight at the gap edge e* = 7T?i/2i: E . A straightforward 
estimate based on the classical density of long trajec- 
tories p(T)dT ~ exp(—T/tu)dT shows that the peak is 
of moderate height p(e*)/ p c \(e*) = 0(1), where p c \(e*) 
is the semiclassical estimate of the DoS. Its width is of 
0(irj/i E ) which reflects an uncertainty in the effective 
trajectory times of 0(£d)- hi a real environment, the 
position of the gap may also be affected by mesoscopic 
fluctuations of system parameters (TBJ. However, such 
effects are beyond the scope of the present paper. 

Summary and discussion — We have solved the quan- 
tum Eilenberger equation to leading order in the small 
parameter tofts- Our solution verifies the existence of 
a gap in the DoS of clean chaotic Andreev billiards. It 
is worthwhile to compare these results to two earlier ap- 
proaches to dealing with the singularities of the classical 
Eilenberger theory: in [7 , Silvestrov et al. argued that 
on bands narrower than a Planck cell, classical dynam- 



ics may be effectively replaced by RMT modeling. In [17 
Vavilov and Larkin, coupled the system to artificial short 
range disorder, fine tuned in strength to mimic quantum 
corrections to classical propagation. This latter proce- 
dure renders the long time dynamics effectively stochas- 
tic, thus preventing the build-up of sharply defined phase 
space structures. Our analysis shows that phenomcno- 
logical input of cither type is not, in fact, necessary. 
The conjunction of classical hyperbolicity and quantum 
uncertainty encoded in the native Eilenberger equation 
automatically regularizes classical singularities at large 
times. This mechanism operates under rather general 
conditions and can be described at moderate theoreti- 
cal efforts. We therefore believe the concepts discussed 
above to be of wider applicability. 

We are grateful for discussions with P. Brouwer. This 
work was supported by SFB/TR 12 of the Deutsche 
Forschungsgemeinschaft and the U.S. Department of En- 
ergy, Office of Science, under Contract No. DE-AC02- 
06CH11357. 



[1] A. Lodder, Y. V. Nazarov, Phys. Rev. B 58, 5783 (1998). 

[2] J. Melsen et al, Europhys. Lett. 35, 7 (1996). 

[3] G. Eilenberger, Z. Phys. 214, 195 (1968); J. Rammer, 
Quantum Field Theory of N on- Equilibrium States (Cam- 
bridge, New York, 2007). 

[4] The distribution of path lengths T in chaotic billiards is 
p{T) oc e~ T/tD 0. 

[5] W. Bauer et al, Phys. Rev. Lett. 65, 2213 (1990). 

[6] L. Wirtz et al, Phy. Rev. B, 56, 7589 (1997). 

[7] P. G. Silvestrov et al, Phys. Rev. Lett. 90, 116801 
(2003). 

[8] R. S. Whitney et al. , Phys. Rev. Lett. 94, 116801 (2005). 
[9] The principal role of the Planck cell in this con- 
text has already been discussed in H. Schomerus and 
C. W. J. Beenakker, Phys. Rev. Lett. 82, 2951 (1999). 

[10] W. H. Zurek, J. P. Paz, Phys. Rev. Lett. 72, 2508 (1994). 

[11] D. J. Gross et al, Adv. Theor. Math. Phys. 4, 893 (2000). 

[12] At its terminal points, G c respects the boundary condi- 
tions 141. This is because at 77 UjSj £) = ±T( u , s ,t)/2 the tra- 
jectorynas either hit the superconductor interface (Fig. [3] 
b,c)), or they have left the vicinity of the reference trajec- 
tory (Fig. [3] a)), meaning that they will hit the interface 
after a classically short time to- 

[13] A detail proof of this property can be found in Appendix 
A of M. Chaichian et al, Phys. Lett. B 666, 199 (2008). 

[14] The (super)linear scaling of V in s implies that as a func- 
tion of u, V(u,s ~ h) < h°. By the same token, the it- 
dependence of the function U(u, s ~ h) ~ h° is classical, 
and the same applies to the mode spectrum. 

[15] This is seen by straightforward comparison of the action 
of the Moyal product's higher order derivatives on the 
envelope function / and its argument function r 9 , resp. 

[16] P. G. Silvestrov, Phys. Rev. Lett. 97, 067004 (2006). 

[17] M. G. Vavilov, A. I. Larkin, Phys. Rev. B 67, 115335 
(2003). 



